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ABSTRACT 


The  viscous  siiock-Iayer  equations  for  three-dimensional 
hypersonic  flows  of  a  perfect  gas  or  equilibrium  air  over  lifting 
bodies  at  high  angles  of  attack  have  been  developed.  For  the  complex 
three-dimensional  reentry  vehicle  geometries  of  interest,  the  re¬ 
sulting  equations  are  written  in  a  nonorthogonal ,  body-oriented 
coordinate  system,  and  the  three  velocity  components  are  defined 
in  the  nonorthogonal  coordinate  directions.  Since  the  viscous 
shock-layer  governing  equations  are  parabolic  in  both  the  stream- 
wise  and  crossflow  directions,  the  equations  are  solved  by  a  highly 
efficient  finite-difference  scheme.  The  principal  advantages  of 
this  technique  are  (i)  the  numerical  method  can  be  used  to  predict 
the  flowfleld  about  arbitrary  geometries  in  both  subsonic  and  super¬ 
sonic  regions,  (ii)  the  solution  is  direct,  and  (ili)  the  effects  of 
Inviscid-vlscous  interactions  are  uniformly  valid  throughout  the 
shock  layer.  Numerical  solutions  have  been  obtained  for  a  1:1. A 
(perfect  gas),  1:2  ellipsoid  in  a  perfect  gas  or  equilibrium  air 
and  the  nose  portion  of  the  shuttle  orbiter  at  zero,  10  and  25-deg 
angles  of  attack.  Comparisons  were  made  with  inviscid  solutions 
and  existing  experimental  data,  and  the  agreement  is  good  for  all 
the  cases. 


INTRODUCTION 

Various  analyses  are  available  to  investigate  the  flowfields  and  wall- 
measurable  quantities  for  hypersonic  flow  past  an  axisymmetric  blunt  body  at 
different  angles  of  attack.  However,  the  effects  of  the  noncircular  body 
geometry  can  be  significant,  and  the  main  purpose  of  this  study  is  to  develop 
a  technique  to  predict  hypersonic  viscous  flows  over  a  noncircular  body  at 
typical  planetary  entry  conditions. 

Recently  a  numerical  method  for  3-D  laminar,  transitional  and/or  turbu¬ 
lent  hypersonic  flows  of  perfect  gas  over  a  blunt  body,  used  for  planetary 
probes,  has  been  Investigated*  at  different  angles  of  attack  by  using  an 
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implicit  finite-difference  viscous  shock-layer  analysis.  However,  that 
analysis  is  valid  for  axisyminetric  bodies  only.  For  the  complex  three- 
dimensional  lifting  reentry  vehicle  geometries  of  interest,  the  resulting 
equations  are  written  in  a  nonorthogonal  body-oriented  coordinate  system, 
and  the  three  velocity  components  are  defined  in  the  nonorthogonal  coordinate 
directions.  This  procedure  is  different  from  writing  the  equations  in  an 
orthogonal  coordinate  system  and  explicitly  performing  a  coordinate  trans¬ 
formation.  Since  the  viscous  shock-layer  equations  are  parabolic  in  both  the 
streamwise  and  crossflow  directions,  the  equations  were  solved  by  a  highly 
efficient  finite-difference  scheme^ ♦^.  The  principle  advantages  of  this 
technique  are  (i)  the  numerical  method  can  be  used  to  predict  the  flowfield 
about  arbitrary  geometries  in  both  subsonic  and  supersonic  regions,  (ii)  the 
solution  is  direct,  and  (ili)  the  effects  of  inviscld-viscous  interactions 
are  Included  within  a  single  set  of  governing  equations  which  are  uniformly 
valid  throughout  the  shock  layer. 

The  basic  formulation  of  the  problem  is  presented  in  Section  2.  Boundary 
conditions  at  the  body  surface  and  immediately  behind  the  shock  are  given  in 
Section  3.  Coordinates  and  method  of  solution  are  presented  in  Sections  4 
and  5,  respectively.  The  results  are  discussed  in  Section  6,  and  the  con¬ 
clusions  are  presented  in  Section  7. 


BASIC  FORMULATION 

The  basic  equations  are  derived  from  the  steady  Navier-Stokes  equations 
in  general  body-oriented  tensor  form  (Figure  1).  One  of  the  coordinates, 

Cl,  is  chosen  in  the  general  axial  direction,  and  another,  C2,  in  a  direction 
normal  to  the  body,  and  the  third,  C3,  around  the  body.  The  coordinate  sys¬ 
tem  requires  orthognality  only  at  the  body  surface.  The  normal  velocity  v 
and  normal  coordinate  C2  are  assumed  to  be  of  order  t  ,  and  all  terms  which  are 
of  higher  order  than  e  are  neglected.  The  methods  of  obtai'.ing  these  equa¬ 
tions  are  discussed  in  detail  elsewhere^*®.  The  nondimensional  form  of  the 
viscous  shock-layer  equations  that  are  applicable  in  the  present  case  can  be 
written  as  follows: 

Continuity  Equation: 


(pug'^)  +  (pvg"^)  +  g|-  (pwg"^)  =  0 


(1) 
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It  should  be  noted  that  the  Prandtl  number  is  assumed  constant  for  a 
perfect  gas  and  is  a  function  of  pressure  and  temperature  for  a  gas  in  chemi¬ 
cal  equilibrium. 

Equation  of  State 

p  =  p(p,  h)  (6) 

For  a  perfect  gas.  Equation  (6)  has  the  analytical  form 

p  =  yp/CCy-i )RT1 


where  y  is  the  ratio  of  specific  heats.  For  the  chemical  equilibrium  as¬ 
sumption,  the  functional  relation  may  be  given  by  a  table  or  an  approximate 
analytical  expression  (curve  fits)  which  are  given  in  detail  in  Reference  4 


In  the  governing  equations  are  Christoffel  symbols  of  the  second 

kind,  and  the  matrix  g..  can  be  obtained  numerically  from  the  grid  generation 
procedure. 


BOUNDARY  CONDITION 


In  order  to  solve  the  above  set  of  governing  equations,  it  is  essential 
to  specify  appropriate  boundary  conditions  at  the  body  surface  and  at  the 
shock.  At  the  body  surface  (wall),  no-sllp  and  no-temperature-jump  conditions 
were  used.  Consequently,  u^  =  Vy  =  w^  =  0,  and  the  wall  temperature  or  heat- 
transfer  rate  was  specified.  The  conditions  immediately  behind  the  shock  were 
obtained  by  using  the  Rankine-Hugonlot  relations.  However,  before  those  equa¬ 
tions  may  be  applied,  the  tangential  and  normal  components  of  the  velocity 
must  be  found.  These  component^  are  obtained  by  the  method  of  Rakich^.  An 
orthogonal  set  of  vectors  with  N  a  unit  vector  normal  to  the  shock  surface  and 
T  a  tangent  vector  is  considered.  The  N-T  plane  is  parallel  to  the  direction 
of  freestream  U^.  Then  the  two-dimensional  Rankine-Hugoniot  relations  are 
used  to  calculate  the  conditions  behind  the  shock.  The  total  velocity  can  be 
written  as 


U  =  T  +  N 


+ 


+ 


(7) 
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where  a,  b,  c,  O  and  6  are  given  by  Raklch^. 
shock  has  been  calculated,  it  is  necessary  to 
body-normal  nonorthogonal  coordinate  system, 
vectors  e^,  e^,  e^  in  terms  of  g^,  in 
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After  the  velocity  behind  the 
rotate  these  components  into  a 
This  is  done  by  expressing  the 
the  ^2*  Cj  directions 


(8) 


Using  Equation  8,  the  shock  velocity  components  in  the  body-normal  non¬ 
orthogonal  coordinate  system  are  determined. 


COORDINATE  GENERATION 

Based  on  the  general  curvilinear  coordinate  governing  equations,  a  body- 
oriented  nonorthogonal  coordinate  system  is  constructed.  This  is  first  done 
on  the  surface  of  the  body,  where  C2  “  0,  and  then  extended  to  the  points 
away  from  the  surface  of  the  body.  The  coordinate  system  requires  ortho¬ 
gonality  only  at  the  body  surface,  and  the  C2  coordinate  is  always  orthogonal 
to  and  ^3.  The  approach  used  is  an  extension  of  that  presented  by  Blottnerf 

and  a  detailed  discussion  of  a  similar  procedure  can  be  found  in  Reference  3. 


METHOD  OF  SOLUTION 

Davis^  presented  an  implicit  finite-difference  method  to  solve  the 
viscous  shock-layer  equations  for  axially  symmetric  flows,  Murray  and  Lewis^ 
extended  the  method  of  solution  to  three-dimensional  high  angle  of  attack 
conditions.  The  present  method  of  solution  is  identical  to  that  of  Murray 
and  Lewis.  Therefore,  only  an  overview  of  the  solution  procedure  is  presented. 

The  solution  begins  on  the  blunt  nose  by  obtaining  an  approximate  stag¬ 
nation  solution  in  the  wind-fixed  coordinate  system.  The  3-D  solution  begins 
on  the  windward  plane  and  marches  around  the  body  obtaining  a  converged  so¬ 
lution  in  each  (j)  sweep.  After  completing  a  sweep  in  the  ())  direction,  the 
procedure  then  steps  downstream  in  and  begins  the  next  <1  sweep.  At  each 
point  the  equations  are  solved  in  the  following  order  (i)  ^i-momentum,  (ii) 
energy,  (ill)  Ci-momentum,  (iv)  integration  of  continuity  for  Y  ,  and  (v) 
coupled  continuity  and  normal  momentum  equations. 


RESULTS  AND  DISCUSSION 

Numerical  solutions  of  the  three-dimensional  nonorthogonal  shock-layer 
equations  were  obtained  for  a  1:1.4,  1:2  ellipsoid  and  the  nose  portion  of 


the  shuttle  orbiter  at  different  angles  of  attack.  It  should  be  mentioned 
here  that  the  wind-fixed  coordinate  system  is  used  in  the  1:2  ellipsoid 
10-deg  angle  of  attack  case.  The  freestream  conditions  for  these  cases  are 
given  in  Table  1.  Comparisons  are  made  for  shock  standoff  distance,  surface 
pressure  distribution  and  heat-transfer  rate  with  Invlscld  results  and  exist¬ 
ing  experimental  data. 

Shock-standoff  distance  as  a  function  of  the  coordinate  along  the  body 
surface  in  different  ({j-planes  is  Illustrated  in  Figures  2  and  3.  The  in- 
viscid  solutions  from  Marconi  and  Yaeger®  are  also  presented  and  are  in  very 
good  agreement  with  the  present  nonorthogonal  viscous  shock-layer  calculations. 
It  is  noted  that  the  shock  standoff  distance  Increases  more  rapidly  between 
((>  =  0  and  30-deg  for  the  1:2  ellipsoid  case  than  for  the  1:1.4  case.  However, 
this  would  be  expected  because  the  body  curvature  in  the  transverse  direction 
is  smaller  for  the  1:2  ellipsoid  case  in  this  region. 

Figures  4  and  5  present  the  surface  pressure  distribution  along  the  body 
at  different  ({/-planes  for  these  two  cases.  Some  Inviscid  pressure  results  are 
also  presented.  It  is  noted  that  the  inviscid  data  are  in  very  good  agree¬ 
ment  with  present  results. 

Figures  6  and  7  show  the  convective  heating  rate  in  several  ({)-planes 
along  the  body  surface.  The  solution  for  Case  1  is  compared  with  the  experi¬ 
mental  data  given  by  Hillsamer  and  Rhudy®.  The  results  show  that  the  present 
prediction  is  in  excellent  agreement  with  the  experimental  data.  Because  of 
the  streamwlse  body  curvature  effects,  the  heating  rate  is  higher  at  (|)  =  90 
and  30-deg  than  ({)  =  0  deg  near  the  stagnation  region  for  the  1:2  ellipsoid 
case.  Except  in  the  stagnation  region,  as  was  expected,  the  convective  heat- 
transfer  rate  decreases  with  increasing  i{). 

Figures  8,  9  and  10  present  the  modified  1:2  ellipsoid  geometry,  surface 
pressure  and  convective  heating  rate  for  an  equilibrium  air  calculation.  As 
shown  in  Figure  9,  the  surface  pressure  is  in  reasonably  good  agreement  with 
the  inviscid  solution.  Since  the  geometry  near  the  stagnation  region  is  very 
flat  and  narrow,  the  convective  heating  rate  decreased  rapidly  between  (|)  =  0 
and  30-deg.  After  ({)  =  30  deg,  a  much  smaller  heating  rate  change  is  noticed 
in  Figure  10. 

Figures  11  to  16  are  the  body  geometry  and  results  of  the  space  shuttle 
orbiter  nose  portion  up  to  seven  nose  radii.  The  cross  section  and  longi¬ 
tudinal  slope  of  the  shuttle  orbiter  are  given  in  Figure  11.  It  is  noted 
that  on  the  windward  side,  the  body  is  relatively  flat  and  the  transverse 
direction  curvature  changes  rapidly  at  about  ({)  =  65-deg.  The  shock  standoff 
distance  in  different  ({/-planes  along  the  body  is  shown  in  Figure  12.  Because 
of  the  flat  windward  side,  the  shock  standoff  distance  decreased  slightly  on 
the  windward  side  and  then  Increased  rapidly  on  the  leeward  side.  The  in¬ 
viscid  results  are  also  presented,  and  the  agreement  is  very  good.  Figure  13 
presents  the  surface  pressure  distribution  around  the  body  at  the  s  =  5.32 
station.  A  very  small  pressure  change  was  detected  on  the  windward  side 
(♦  —  50°).  However,  a  rapidly  decreasing  surface  pressure  was  noticed  in  the 
region  50  <  (^<90  due  to  the  transverse  curvature  effect.  On  the  leeward 

side  (4/  >  90°),  the  pressure  tends  to  stabilize  and  the  change  is  small  once 
again.  Figure  14  shows  that  the  same  trend  is  followed  for  all  stations 
along  the  body  surface.  Finally,  the  temperature  profiles  normal  to  the  body 


surface  at  s  =  0  (stagnation  point)  and  s  =  5.325  and  the  convective  heating 
rate  along  the  body  surface  are  plotted  in  Figures  15  and  16.  Since  no  other 
dat^>  are  available  for  this  case,  only  the  present  results  are  presented  here 
For  the  same  reason  discussed  earlier,  the  temperature  and  heating  rate  jump 
in  the  region  60  <  <|)  <  90  was  also  noticed. 


CONCLUSIONS 

The  main  objective  of  this  study  is  to  investigate  the  influence  of  the 
noncircular  body  on  the  shock-layer  flow  phenomena  within  reasonable  com¬ 
puting  times.  A  general  nonorthogonal  shock-layer  analysis  is  used.  The 
present  surface  pressure,  shock  standoff  distance  and  convective  heating  rate 
predictions  are  in  good  agreement  with  inviscid  solutions  and  existing  ex¬ 
perimental  data.  The  viscous  shock-layer  technique  requires  much  less  computing 
time  than  the  time-dependent  method*®  (Table  1).  In  the  future  we  plan  to 
Include  the  transitional  and/or  turbulent  flow  in  the  technique. 
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LIST  OF  SYMBOLS 

specific  heat  at  constant  pressure 

unit  vector  of  a  general  orthogonal  coordinate  system 

vectors  in  (streamwise) ,  ^2  (normal)  and  C3  (circumferential) 
directions 

.  .  _  .  dr  3r  .  2  9d>  96 

metrics  gy  j  +  j  +  r 
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^ij^ 
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22  ^13 


* ,  *2 

static  enthalpy,  h  /U 


shock-normal  vector 
pressure 
Prandtl  number 
body  nose  radius 


temperature 

e-  *2 ,  * 

reference  temperature,  U  /C 

CO  p 


shock-tangent  vector 


i. 


streamwise,  normal  and  crossflow  velocity  components  nondimen- 
sionalized  by  the  freestream  velocity  U* 
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"^ABLE  1. 

Test  Case  Conditions 

Case 

1 

2 

3 

4 

Body  Geometry 

1:1. A 

(ellipsoid) 

1:2 

(ellipsoid) 

1:2 

(modified 

ellipsoid) 

Shuttle 

o 

81.5 

203.5 

540.0 

94,87 

Poo»  slug/ft^ 

4.118E-6 

2 . 3A0E-A 

7.99E-4 

1.0528E-3 

P„.  psia 

O.OAl 

0.567 

5.145 

0.0119 

T  (°R) 
w 

A70 

366 

3139 

540 

Re  (1/ft) 

00 

2.38E-f6 

4.993E46 

3.3058E47 

5.5696E45 

a,  deg 

0° 

10° 

0° 

25° 

Gas  Model 

perfect  gas 

perfect  gas 

equil ibrium 
air 

perfect  gas 

Computational 
Time* (Min) 

IBM  370/3032 

6 

12 

16 

12 

*It  takes  approximately  two  hours  CDC  7600  time  for  a  sphere  nose  when  using 
time  dependent  method  of  Reference  10. 
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Figure  1.  Coordinate  System 


Figure  5.  Surface  Pressure  for  1:2  Ellipsoid 


Figure  6.  Convective  Heat  Transfer  to  1:43  Ellipsoid. 
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Figure  9.  Surface  Pressure  for  1:2  Ellipsoid 


Flgurp  10.  Convective  ile.nt  Transfer  to  1:2  Ellipsoid. 
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Figure  15.  Space  Shuttle  Temperature  Profile  at  Stagnation 
Point  and  s  =  5.32 
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